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The non-equilibrium dynamics of mixing, oscillations and equilibration is studied in a field theory 
of flavored neutral mesons that effectively models two flavors of mixed neutrinos, in interaction with 
other mesons that represent a thermal bath of hadrons or quarks and charged leptons. This model 
describes the general features of neutrino mixing and relaxation via charged currents in a medium. 
The reduced density matrix and the non-equilibrium effective action that describes the propagation 
of neutrinos is obtained by integrating out the bath degrees of freedom. We obtain the dispersion 
relations, mixing angles and relaxation rates of "neutrino" quasiparticles. The dispersion relations 
and mixing angles are of the same form as those of neutrinos in the medium, and the relaxation 
rates are given by Ti(fc) = F ee (k) cos 2 9 m (k) + r MM (fc) sin 2 6 m (k) ; r 2 (fc) = r MM (fc) cos 2 8 m (k) + 
r ee (fc) sin 2 9 m (k) where T aa (k) are the relaxation rates of the flavor fields in absence of mixing, 
and 9 m (k) is the mixing angle in the medium. A Weisskopf-Wigner approximation that describes 
the asymptotic time evolution in terms of a non-hermitian Hamiltonian is derived. At long time 
>> "neutrinos" equilibrate with the bath. The equilibrium density matrix is nearly diagonal 
in the basis of eigenstates of an effective Hamiltonian that includes self-energy corrections in the 
medium. The equilibration of "sterile neutrinos" via active-sterile mixing is discussed. 

PACS numbers: 13.15.+g,12.15.-y,f l.fO.Wx 



I. INTRODUCTION 



Neutrinos are the central link between particle and nuclear physics, astrophysics and cosmology [E HI H> & Q and the 
experimental confirmation of neutrino mixing and oscillations provide a first evidence for physics beyond the Standard 
Model. Neutrino mixing provides an explanation for the solar neutrino problem^, 0, &H, plays a fundamental role in 



15, [16|, [17| and in cosmology 1 18] : in big bang nucleosynthesis 
231 ] . structure forrnation[24l |25|. 1261 1271 128|| and possible dark 



the physics of core-collapse supernovae [1 0. 1 1 ll. 1 1 2L 1 1 c 
(BBN)[19||, baryogenesis through leptogenesis [20. |21|. 
matter candidate Hi, [27|, [H[ . 

The non- equilibrium dynamics of neutrino mixing, oscillations and equilibration is of fundamental importance in all 
of these settings. Neutrinos are produced as "flavor eigenstates" in weak interaction vertices, but propagate as a linear 
superposition of mass eigenstates. This is the origin of neutrino oscillations. Weak interaction collisional processes 
are diagonal in flavor leading to a competition between production, relaxation and propagation which results in a 
complex and rich dynamics. 

Beginning with pioneering work on neutrino mixing in media [29l. [30l l3ll l32l . |33| . the study of the dynamical 
evolution has been typically cast in terms of single particle "flavor states" or matrix of densities that involve either 
a non-relativistic treatment of neutrinos or consider flavor neutrinos as massless. The main result that follows from 
these studies is a simplified set of Bloch equations with a semi-phenomenological damping factor (for a thorough 
review see[18|). 

Most of these approaches involve in some form the concept of distribution functions for "flavor states" , presumably 
these are obtained as expectation values of Fock number operators associated with flavor states. However, there are 
several conceptual difficulties associated with flavor Fock states, still being debated [3~S, l35l l36l l37l. l38l l39l. |40L l4l| ■ 

The importance of neutrino mixing and oscillations, relaxation and equilibration in all of these timely aspects of 
cosmology and astroparticle physics warrant a deeper scrutiny of the non-equilibrium phenomena firmly based on 
quantum field theory. 

The goals of this article: Our ultimate goal is to study the non-equilibrium dynamics of oscillation, relaxation 
and equilibration directly in the quantum field theory of weak interactions bypassing the ambiguities associated 
with the definition of flavor Fock states. We seek to understand the nature of the equilibrium state: the free field 
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Hamiltonian is diagonal in the mass basis, but the interactions are diagonal in the flavor basis, however, equilibration 
requires interactions, hence there is a competition between mass and flavor basis, which leads to the question of which 
is the basis in which the equilibrium density matrix is diagonal. Another goal is to obtain the dispersion relations 
and the relaxation rates of the correct quasiparticle excitations in the medium. 

In this article we make progress towards these goals by studying a simpler model of two "flavored" mesons repre- 
senting the electron and muon neutrinos that mix via an off-diagonal mass matrix and interact with other mesons 
which represent either hadrons (neutrons and protons) or quarks and charged leptons via an interaction vertex that 
models the charged current weak interaction. The meson fields that model hadrons (or quarks) and charged leptons 
are taken as a bath in thermal equilibrium. In the standard model the assumption that hadrons (or quarks) and 
charged leptons can be considered as a bath in thermal equilibrium is warranted by the fact that their strong and 
electromagnetic interactions guarantee faster equilibration rates than those of neutrinos. 

This model bears the most relevant characteristics of the standard model Lagrangian augmented by an off diagonal 
neutrino mass matrix and will be seen to yield a remarkably faithful description of oscillation and relaxational 
dynamics in a thermal medium at high temperature. It effectively describes the thermalization dynamics of neutrinos 
in a medium at high temperature such as the early Universe for T > 3 MeVpjL [Til f]~o(] . 

Furthermore, Dolgov et. al.[13] argue that the spinor nature of the neutrinos is not relevant to describe the dynamics 
of mixing at high energies, thus we expect that this model captures the relevant dynamics. 

An exception is the case of neutrinos in supernovae, a situation in which neutrino degeneracy, hence Pauli blocking, 
becomes important and requires a full treatment of the fcrmionic aspects of neutrinos. Certainly the quantitative 
aspects such as relaxation rates must necessarily depend on the fermionic nature. However, we expect that a bosonic 
model will capture, or at minimum provide a guiding example, of the most general aspects of the non-equilibrium 
dynamics. The results found in our study lend support to this expectation. 

While meson mixing has been studied previously [43j, mainly motivated by mixing in the neutral kaon and pseu- 
doscalar 77, rf systems, our focus is different in that we study the real time dynamics of oscillation, relaxation and 
equilibration in a thermal medium at high temperature including radiative corrections with a long view towards un- 
derstanding general aspects that apply to neutrino physics in the high temperature environment of the e arly Universe. 

While neutrino equilibration in the early Universe for T > 3MeV prior to BBN is undisputable[3l. [l8l 1190, the mam 
questions that we address in this article arc whether the equilibrium density matrix is diagonal in the flavor or mass 
basis and the relation between the relaxation rates of the propagating modes in the medium. 

The strategy: the meson fields that model flavor neutrinos are treated as the "system" while those that describe 
hadrons (or quarks) and charged leptons, as the "bath" in thermal equilibrium. An initial density matrix is evolved 
in time and the "bath" fields are integrated out up to second order in the coupling to the system, yielding a "reduced 
density matrix" which describes the dynamics of correlation functions solel y of system fields (neutrinos). This program 
pioneered by Feynman and Vernon[44| for coupled oscillators (see also (45l |46| ) is carried out in the interacting theory 
by implementing the closed-time path-integral re pres entation of a time evolved density matrix |48|. This method 
yields the real time non- equilibrium effective action 47] including the self-energy which yields the "index of refraction" 
correction to the mixing angles and dispersion relations [4!| in the medium and the decay and relaxation rates of the 
quasiparticle excitations. The non-equilibrium effective action thus obtained yields the time evolution of correlation 
and distribution functions and expectation values in the reduced density matrix [47j. The approach to equilibrium is 
determined by the long time behavior of the two point correlation function and its equal time limit, the one-body 
density matrix. The most general aspects of the dynamics of mixing and equilibration are completely determined by 
the spectral properties of the correlators of the bath degrees of freedom in equilibrium. 

Brief summary of results: 

• We discuss the ambiguities in the definition of flavor Fock operators, states and distribution functions. 

• The non-equilibrium effective action is obtained up to second order in the coupling G ~ Gf between the "system" 
(neutrinos) and the bath (hadrons, quarks and charged leptons) in equilibrium. It includes the one-loop matter 
potential contribution (0(G)) and the two-loop (C(G 2 )) retarded self-energy. The "index of refraction" (49j 
is determined by the matter potential and the real part of the space-time Fourier transform of the retarded 
self-energy. The relaxation rates of the quasiparticle excitations are determined by its imaginary part. The 
non-equilibrium effective action leads to Langevin-like equations of motion for the fields with a noise term 
determined by the correlations of the bath, it features a Gaussian probability distribution but is colored. The 
noise correlators and the self-energy fulfill a generalized fluctuation-dissipation relation. 

• We obtain expressions for the dispersion relations and mixing angles in medium which are of the same form as 
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in the case for neutrinos. The relaxation rates for the two types of quasiparticles are given by 

ri(fc) = T ee (k)cos 2 e m (k) + r^(k)sm 2 e m (k) (1.1) 
r 2 (fc) = r MM (fc) cos 2 e m {k) + v ee {k) sin 2 e m {k) (1.2) 

where T aa (k) are the relaxation rates of the flavor fields in absence of mixing, and 9 m {k) is the mixing angle in 
the medium. 

• A Weisskopf-Wigner description of the long time dynamics in terms of an effective non-hermitian Hamiltonian 
is obtained. Although this effective description accurately captures the asymptotic long time dynamics of the 
expectation value of the fields in weak coupling, it does not describe the process of equilibration. 

• For long time >> T^J, the two point correlation function of fields becomes time translational invariant reflect- 
ing the approach to equilibrium. The one-body density matrix reaches its equilibrium form at long time, in 
perturbation theory it is nearly diagonal in the basis of eigenstates of an effective Hamiltonian that includes self- 
energy corrections in the medium, with perturbatively small off-diagonal corrections in this basis. The diagonal 
components are determined by the distribution function of eigenstates of this in- medium effective Hamiltonian. 

• These results apply to the case of sterile neutrinos with modifications to the dispersion relations and relaxation 
rates arising from simple "sterility" conditions. "Sterile neutrinos" equilibrate with the bath as a consequence 
of active-sterile mixing [U [H, [5(| ■ 

In section |TIJ we introduce the model and discuss the ambiguities in defining flavor Fock operators, states and 
distribution functions. In section (|IIIj) we obtain the reduced density matrix, the non-equilibrium effective action and 
the Langevin-like equations of motion for the expectation value of the fields. In section (|IV[) we provide the general 
solution of the Langevin equation. In section (|V| we obtain the dispersion relations, mixing angles and decay rates 
of quasiparticle modes in the medium. In this section an effective Weisskopf-Wigner description of the long time 
dynamics is derived. In section (|VIj) we study the approach to equilibrium in terms of the one-body density matrix. 
In this section we discuss the consequences for "sterile neutrinos" . Section (jVIIj) summarizes our conclusions. 



II. THE MODEL 



We consider a model of mesons with two flavors e , fj, in interaction with a "charged current" denoted W and 
a "flavor lepton" Xa modeling the charged current interactions in the electroweak (EW) model. In terms of field 
doublets 



the Lagrangian density is 



f 



£ = ~ {d^ T 8^ - $ T M 2 $} + C [W,x] + GW^ T -X + G&xl + G<gxl (2.2) 



where the mass matrix is given by 



M 2 M 2 
M 2 M 2 



(2.3) 



where £o[Wi x) is the free field Lagrangian density for W, x which need not be specified. The mesons <f> eyfl play the role 
of the flavored neutrinos, Xe,^ the role of the charged leptons and W a charged current, for example the proton-neutron 
current Jry p (l — gA"f5)n or a similar quark current. The coupling G plays the role of Gf- As it will be seen below, 
we do not need to specify the precise form, only the spectral properties of the correlation function of this current are 
necessary. 

Passing from the flavor to the mass basis for the fields cf> e ^ by an orthogonal transformation <!> = U(0) ip 

U ® I £ I 5 m = ( ^ Zl ) (2.4) 
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where the orthogonal matrix U{6) diagonalizes the mass matrix M 2 , namely 



U- l {6) M 2 U(d) 



Ml 
Mf 



(2.5) 



In the flavor basis M can be written as follows 



where we introduced 



I 2 = 17 2 1 + — 2 ^ ~~ cos 26 sin 29 



sin 29 cos 26* 



(2.6) 



m 2 = \ i M l + M I) ; sm 2 = M l - M l 



(2.7) 



Mass and flavor states: 



It is convenient to take the spatial Fourier transform of the fields <j) a ; ipi and their canonical momenta tt 
ipi with a = e, ji and i = 1, 2 and write (at t=0), 



0a! 



4>a(x) 



1 



Ak-x 



(2.8) 



in these expressions we have denoted the spatial Fourier transforms with the same name to avoid cluttering of notation 
but it is clear from the argument which variable is used. The free field Fock states associated with mass eigenstates 
are obtained by writing the fields which define the mass basis ipi in terms of creation and annihilation operators, 



,-k 



with 



y/2ui(k) 



(k) = Jk 2 + M? ; i= 1,2 



i, — k 



(2.9) 



(2.10) 



The annihilation (a i g) and creation (at -) operators obey the usual canonical commutation relations, and the free 
Hamiltonian in the mass basis is the usual sum of independent harmonic oscillators with frequencies u>i(k). One can, 
in principle, define annihilation and creation operators associated with the flavor fields a a g, a - respectively in a 
similar manner 



y/2Sl a (k) 



y/2Sl a (k) 



(2-11) 



with the annihilation (o Q g) and creation (a -) operators obeying the usual canonical commutation relations. How- 
ever, unlike the case for the mass eigenstates, the frequencies £l a (k) are arbitrary. Any choice of these frequencies 
furnishes a different Fock representation, therefore there is an intrinsic ambiguity in defining Fock creation and an- 
nihilation operators for the flavor fields since these do not have a definite mass. In references^, [H, Hi| a particular 
assignment of masses has been made, but any other is equally suitable. The orthogonal transformation between the 
flavor and mass fields eqn. (]2 .4[) . leads to the following relations between the flavor and mass Fock operators, 
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V£ = cosf 



a 1 s A e4 (fc) +aj _-B e ,i(fc) 



sin ( 



sin 



a 2 ^A e . 2 (k) + a* _-B e!2 (A;) 



(2.12) 



(2.13) 



where A a) i,B a) j are the generalized Bogoliubov cocfhcicnts 



B, 



(2.14) 



These coefficients obey the condition 



(Al ti - Bl ti ) = 1 



(2.15) 



which guarantees that the transformation between mass and flavor Fock operators is formally unitary and both 
sets of operators obey the canonical commutation relations for any choice of the frequencies £l a (k). Neglecting the 
interactions, the ground state |0 > of the Hamiltonian is the vacuum annihilated by the Fock annihilation operators 
of the mass basis, 



a. AO >= for all i = 1,2, k. 



(2.16) 



In particular the number of flavor Fock quanta in the non-interacting ground state, which is the vacuum of mass 
eigenstates is 



<°\ a U a e,k\°> = C0S2 



<°\ a lt a a.k\°> = C0S2 



fl e (k) — u>i(k) 



AQ e (k)uji{k) 
Qftik) - w 2 (fc) 



40 (U (fc)w 2 (A;) 



i 2 



sin 2 9 



An e {k)u 2 {k) 



4n p (*;)wi(fc) 



(2.17) 



(2.18) 



namely the non-interacting ground state (the vacuum of mass eigenstates) is a condensate of "flavor" states 36, 38, 39] 
with an average number of "flavored particles" that depends on the arbitrary frequencies £l a (k). Therefore these 
"flavor occupation numbers" or "flavor distribution functions" are not suitable quantities to study equilibration. 
Assuming that £l a (k) — > k when fc — > oo, in the high energy limit A — > 1 ; B — > and in this high energy limit 



sin 9 a 2 g 



cos 9a 2 g 



sin^aj^ g 



(2.19) 



Therefore, under the assumption that the arbitrary frequencies f2 Q (fc) — > k in the high energy limit, there is an 
approximate identification between Fock states in the mass and flavor basis in this limit. However, such identification 
is only approximate and only available in the asymptotic regime of large momentum, but becomes ambiguous for 
arbitrary momenta. In summary the definition of flavor Fock states is ambiguous, the ambiguity may only be 
approximately resolved in the very high energy limit, but it is clear that there is no unique definition of a flavor 
distribution function which is valid for all values of momentum k and that can serve as a definite yardstick to study 
equilibration. Even the non-interacting ground state features an arbitrary number of flavor Fock quanta depending on 
the arbitrary choice of the frequencies fi, a (k) in the definition of the flavor Fock operators. This is not a consequence 
of the meson model but a general feature in the case of mixed fields with similar ambiguities in the spinor case 4l| . 

We emphasize that while the flavor Fock operators are ambiguous and not uniquely defined, there is no ambiguity 
in the flavor fields <p a which are related to the mass fields (pi via the unitary transformation (|2.4[) . While there is 
no unambiguous definition of the flavor number operator or distribution function, there is an unambiguous number 
operator for the Fock quanta in the mass basis Ni(k) = a ~a i g, whose expectation value is the distribution function 
for mass Fock states. 
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III. REDUCED DENSITY MATRIX AND NON-EQUILIBRIUM EFFECTIVE ACTION 

Our goal is to study the equilibration of neutrinos with a bath of hadrons or quarks and charged leptons in thermal 
equilibrium at high temperature. This setting describes the thermalization of neutrinos in the early Universe prior to 
BBN, for temperatures T > 3 MeV|l, M, El- 

We focus on the dynamics of the "system fields" , either the flavor fields <f> a or alternatively the mass fields tp i . The 
strategy is to consider the time evolved full density matrix and trace over the bath degrees of freedom x, W . It is 
convenient to write the Lagrangian density (|2.2|) as 

C[4> a ,Xa, W] = Co[4>] + Co\W,x\ + G<p a O a + G4>lxl (3.1) 

with an implicit sum over the flavor label a = e, ji, where 

O a = Xa W. (3.2) 

£()[■■■] are the free Lagrangian densities for the fields <fi a ,Xa,W respectively. The fields 4> a are considered as the 
"system" and the fields Xa, W are treated as a bath in thermal equilibrium at a temperature T = 1//3. We consider 
a factorized initial density matrix at a time t% — of the form 

p(Q) = fo{V)®e~ f3Hobc ' W] (3-3) 

where Hq\\, W] is Hamiltonian for the fields x^ ' ■ Although this factorized form of the initial density matrix leads 
to initial transient dynamics, we are interested in the long time dynamics, in particular in the long time limit. The 
bath fields Xa , W will be "integrated out" yielding a reduced density matrix for the fields (f> a in terms of an effective 
real-time functional, known as the influence functional 44] in the theory of quantum brownian motion. The reduced 
density matrix can be represented by a path integral in terms of the non-equilibrium effective action that includes 
the influence functional. This method has been used extensively to study quantum brownian motion [44L l45j . and 
quantum kinetics [46l |47|. 

In the flavor field basis the matrix elements of p<p(0) are given by 

(0alMO)|$}=P*o(0a;$») ( 3 - 4 ) 

or alternatively in the mass field basis 

{<Pi\pip(0)\<Pj) = P<p;o(fi; <Pj) ■ (3.5) 
The time evolution of the initial density matrix is given by 

P (t f ) = e-^^f-^piuy 11 ^?-^ , (3.6) 

where the total Hamiltonian H is 

H = H [cl>]+H \ X ,W}+H I [ ( t>,x,W}. (3.7) 

The calculation of correlation functions is facilitated by introducing currents coupled to the different fields. Further- 
more since each time evolution operator in eqn. (|3.6p will be represented as a path integral, we introduce different 
sources for forward and backward time evolution operators, referred to as J + , J~ respectively. The forward and 
backward time evolution operators in presence of sources are U(tf, U; J + ), U~ 1 {tf , U, J~) respectively. 

We will only study correlation functions of the "system" fields <j) (or ip in the mass basis), therefore we carry out 
the trace over the x and W degrees of freedom. Since the currents allow us to obtain the correlation functions 
for any arbitrary time by simple variational derivatives with respect to these sources, we take tf — > oo without loss of 
generality. The non-equilibrium generating functional is given bv(46l. [4?| 



Z\j + ,T] =TrU(oc,t i ;J+)p(t i )U-\oc,t i ,J-). (3.8) 

Where stand collectively for all the sources coupled to different fields. Functional derivatives with respect to the 
sources J + generate the time ordered correlation functions, those with respect to J~ generate the anti-time ordered 
correlation functions and mixed functional derivatives with respect to J + , J~ generate mixed correlation functions. 
Each one of the time evolution operators in the generating functional (|3.8p can be written in terms of a path integral: 
the time evolution operator U(oo,ti; J + ) involves a path integral forward in time from U to t — oo in presence of 
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sources J + , while the inverse time evolution operator U~ 1 (oo,ti, J") involves a path integral backwards in time from 
t = oo back to t$ in presence of sources J~. Finally the equilibrium density matrix for the bath e-P H a\xW] can ^ e 
written as a path integral along imaginary time with sources J' 3 . Therefore the path integral form of the generating 
functional (|3.8[) is given by 



z\j + ,r 



(3.9) 



with the boundary conditions <i> + (a;, U) = ; $~(x,ti) = &i(x). The trace over the bath fields %, W is performed 

with the usual periodic boundary conditions in Euclidean time. 
The non-equilibrium action is given by 



S[$ ± )X ± ;J£;J±] J%] = 



dtd 3 a 



£ (<t> + ) + J£<t> + -C (<l>-)- J7<t>~ 



l x[c [ X , W] + J xX + JwW + G<p a O a + G^lxl) 



(3.10) 



where C describes the following contour in the complex time plane: along the forward branch (ij,+oo) the fields 
and sources are <fr + , x + , «/+, along the backward branch (oo, U) the fields and sources are <&~, x~, JZ and along the 
Euclidean branch (U,ti — ifj) the fields and sources are $ = 0;%^, Along the Euclidean branch the interaction 
term vanishes since the initial density matrix for the field x is assumed to be that of thermal equilibrium. This contour 
is depicted in fig. flU 




ti-i/3 

FIG. 1: Contour in time for the non-equilibrium path integral representation. 



The trace over the degrees of freedom of the x field with the initial equilibrium density matrix, entail periodic 
boundary conditions for x> W along the contour C. However, the boundary conditions on the path integrals for the 
field $ are given by 



$+(x,t = oo) = $ _ (x,t = oo) (3.11) 

and 

<Z> + (x,t = ti) = $i(x) ; $-(x,t = ti) =%{x) (3.12) 

The reason for the different path integrations is that whereas the x an d W fields are traced over with an initial 
thermal density matrix, the initial density matrix for the $ field will be specified later as part of the initial value 
problem. The path integral over Xi W leads to the influence functional for $ ± [i3|. 

Because we are not interested in the correlation functions of the bath fields but only those of the "system" fields, 
we set the external c-number currents J x = 0; Jw = 0. Insofar as the bath fields are concerned, the system fields $ 
act as an external c-number source, and tracing over the bath fields leads to 

J V^V^VxfiVWPe 10 **^™* ** *^*} = ^/c^W^^-^W. (3.13) 
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The expectation value in the right hand side of eqn. (|3 . 1 3|) is in the equilibrium free field density matrix of the 
fields Xj W. The path integral can be carried out in perturbation theory and the result exponentiated to yield the 
effective action as follows 

d 4 x [ d 4 x'(t> a (x)(t>fs[x')(o a (x)Of3(x')) + 0(G 3 ) (3.14) 
c Jc x ' 



(iG) 



This is the usual expansion of the exponential of the connected correlation functions, therefore this series is identified 
with 

e iGj c d i x^ a o a +4>lx 2 a \ _ e » Li f {</>+, <!>-} / 3 15 \ 

/ 

where Lif[(f> + , <f>~] is the influence functional[4:4\, and (• • • )o stand for expectation values in the bath in equilibrium. 
For (x.a(x)W(x)^ — the influence functional is given by 



L 



■ f [4> + ,<t>-] = G JcPx&ix) (xl(x)) o + i^- J d*x Jd*^M^)M^){O a (x)O p {a/)) + 0(G 3 ) . (3.16) 

In the above result we have neglected second order contributions of the form G 2 </> 4 . These non- linear contributions 
give rise to interactions between the quasiparticles and will be neglected in this article. Here we are primarily concerned 
with establishing the general properties of the quasiparticles and their equilibration with the bath and not with their 
mutual interaction. As in the case of mixed neutrinos, the inclusion of a "neutrino" background may lead to the 
phenomenon of non-linear synchronization}^!! , but the study of this phenomenon is beyond the realm of this 

article. 

We focus solely on the non-equilibrium effective action up to quadratic order in the "neutrino fields" , from which 
we extract the dispersion relations, relaxation rates and the approach to equilibrium with the bath of the quasiparticle 
modes in the medium. 

The integrals along the contour C stand for the following expressions: 

Gj^x^ a (x)(^ a (x)) Q =Vaa j d*x dt[4>l + (x)-<t>l-(x)) (3.17) 

where V aa are the "matter potentials" which are independent of position under the assumption of translational 
invariance, and time independent under the assumption that the bath is in equilibrium, and 



d A x j^x'<t> a (x)<t> (3 (x')(O a (x)O (} (x')) = J d 3 x J dt J d 3 x' J dt'[ct>i(x)ct>+(x')(OZ(x)O+(x')) 

+ cj>Z(x)cj>e(x>)(O-(x)Op(x>)) - 4>+(x)<t>p(x')(O+(x)Op(x>)) 

- 0-(x)4(a/)<O-( a! )O+(a/))o] (3.18) 

Since the expectation values above are computed in a thermal equilibrium translational invariant density matrix, 
it is convenient to introduce the spatial Fourier transform of the composite operator O in a spatial volume V as 

O a % {t) = ^=J d 3 xe ins O a (x,t) (3.19) 
in terms of which we obtain following the correlation functions 

i°:,k^ t-^ = Tt °P-^ r~ PH ° [x > W] O ark (t) = G> p (k; t - = G-+(k; t - t>) (3.20) 

(Ojg(t)f-_g(f)) = TrO Qjfc -(t) e -" ff * <Vg(f ) = G< p {k;t - f) = G^(k;t- f) = ^+(fc;t' - t) (3.21) 

(°tk {t) °i-k {t 'V = Oe(* - + G<p(k;t- t')Q(t> -t) = Q++(k;t- t>) (3.22) 

( Z,H®°l-sW = ^ p (k;t- t')G(t' -t) + G<p(k;t- t>)0(t - t') = G^(k;t- t') (3.23) 
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The time evolution of the operators is determined by the Heisenberg picture of l?o[x, W]- Because the density 
matrix for the bath is in equilibrium, the correlation functions above are solely functions of the time difference as 
made explicit in the expressions above. These correlation functions are not independent, but obey 



0#"(*5 *, + G~p(k; t, t') ~ g-+(k; t, t') - g+-(k ; t, t') = 
The correlation function up to lowest order in the coupling G is given by 

d 3 p 



(3.24) 



(3.25) 



where the expectation value is in the free field equilibrium density matrix of the respective fields. This correlation 
function is diagonal in the flavor basis and this entails that all the Green's functions (|3. 20113.23)1 are diagonal in the 
flavor basis. 

The non-equilibrium effective action yields the time evolution of the reduced density matrix, it is given by 



L 



dtd 3 x [£ o (0 + )-A>(<n] +-M<£ + ,<n 



(3.26) 



where we have set the sources for the fields (jr^ to zero. 

In what follows we take ti — without loss of generality since (i) for t > ti the total Hamiltonian is time independent 
and the correlations will be solely functions of t — ti, and (ii) we will be ultimately interested in the limit t ^> t i when 
all transient phenomena has relaxed. Adapting the methods presented in ref. [471 ] to account for the matrix structure 
of the effective action, introducing the spatial Fourier transform of the fields (fr^ defined as in eqn. (|3.19|) and the 
matrix of the matter potentials 



V 



V ee 

K, 



(3.27) 



we find 



>L 



effl 



E 5 



dt 



G 2 , , 



dt' 



!! jo 

+- 



(3.28) 



The "matter potentials" V aa play the role of the index of refraction correction to the dispersion relations [49j and 
is of first order in the coupling G whereas the contributions that involve Q are of order G 2 . As it will become clear 
below, it is more convenient to introduce the Wigner center of mass and relative variables 



1 



* a @,t) = -(4>+(2,t) + <f>-(3,t)) ; R a (S,t) = (4£@,t) - <f>-(2,t)) 
and the Wigner transform of the initial density matrix for the fields 



(3.29) 



with the inverse transform 



(3.30) 



2 2 J 



R° 



(3.31) 
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The boundary conditions on the <\> path integral given by (|3.12p translate into the following boundary conditions 
on the center of mass and relative variables 



tf a (2,t = 0) = tf° ; R a (x,t = 0) = R° a (3.32) 
furthermore, the boundary condition (|3.11|) yields the following boundary condition for the relative field 

R a (x,t = oo) = 0. (3.33) 

This observation will be important in the steps that follow. 

The same description applies to the fields in the mass basis. We will treat both cases on equal footing with the 
notational difference that Greek labels a, (3 refer to the flavor and Latin indices i,j refer to the mass basis. 

In terms of the spatial Fourier transforms of the center of mass and relative variables (|3.29p introduced above, 
integrating by parts and accounting for the boundary conditions (|3 . 32[) . the non-equilibrium effective action (|3.28[) 
becomes: 



dt 



£ dt J™ dt' £ { t - + ^(k; t ~ O^gCo} 



+ i / d 3 xR° a (x)4> a (x,t = 0) 



(3.34) 

where the last term arises after the integration by parts in time, using the boundary conditions (|3.32p and (|3.33p . 
The kernels in the above effective Lagrangian are given by (see eqns. (|3. 20113. 23"|) ) 

JC af3 (k;t-t') = ^-\G>p(k;t-t') + g< (k;t-t') 



i^Jk-t-t') 



G 



(3.35) 

g>p(k; t-t')- g<p(k; t - t')J 0(i - = iZ a0 (k; t - t')Q(t - t') (3.36) 
The term quadratic in the relative variable R can be written in terms of a stochastic noise as 
exp{ - \ J dt J dt'R a _j:(t)lC a0 (k;t-t')R (t')} = J exp{-^jdtj dt' ^ £ (i)/C^(fc; t - 0^,_fe(*') 

(3.37) 



(3.38) 



J *£a,-S(*)tf«,iE(*)} 
The non-equilibrium generating functional can now be written in the following form 

Z = J J DU° J V^VRVi VV(^r ; n )D J R e^ d3a:iJo '«^ ) ( n °(*)-*«(^ i = ))p[^ 



exp 



dtR a dt) 



+ ( fc2 <W + M '/3 + Va^gW + / dt' K P {k- 1 - 0*^(0 - £ ati {t) 



(3.39) 



The functional integral over R° can now be done, resulting in a functional delta function, that fixes the boundary 
condition ^f a (x,t = 0) = 11° (x). Finally the path integral over the relative variable can be performed, leading to a 
functional delta function and the final form of the generating functional given by 



Z = J D^ DU W(^ a ;U a )V^V^V[^} x 

S K,g(*) + ( fc2 S °P + M lp + Y ^^(t) + I dt' £ Q/3 (fc; * - t')* % {t') - t a % {t) 



with the boundary conditions on the path integral on given by 

* a (2,i = 0) ; i> a (x,t = 0) = Il o a {x) 



(3.40) 



(3.41) 
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where we have used the definition of S?g(/c; t — t') in terms of S Q ^(A:; t — t') given in equation (|3.36[) . 

The meaning of the above generating functional is the following: to obtain correlation functions of the center of 
mass Wigner variable ^ we must first find the solution of the classical stochastic Langevin equation of motion 

r t 

'„2 x . i H>*2 



+ ( fc S aP + M '/5 + Va/jJ^.gW + / dt ' X a p(k; t - H)V pjtf) = ^(t) 

w 

= o) = *a,s ; * Q ,^ = °) = n L- ( 3 - 42 ) 

for arbitrary noise term £ and then average the products of iff over the stochastic noise with the Gaussian probability 
distribution V[£] given by (|3.39[) . and finally average over the initial configurations *$f°(x);Il (x) weighted by the 
Wigner function VV(\^°, II ), which plays the role of an initial phase space distribution function. 

Calling the solution of (I3.42|) iff ^IT), the two point correlation function, for example, is given by 



_ f v[€\P[£] f dv° f dip w(w°;n°)w ai g(t;C;^ ;n°)^_g(t^;^ ;n ) 



In computing the averages and using the functional delta function to constrain the configurations of \& to the solutions 
of the Langevin equation, there is the Jacobian of the operator (d 2 /dt 2 + k 2 )S a p + M 2 + Y + J dt'Y, ret (k; t — t') which 
however, is independent of the field and the noise and cancels between numerator and denominator in the averages. 
There are two different averages: 

• The average over the stochastic noise term, which up to this order is Gaussian. We denote the average of a 
functional over the noise with the probability distribution function P[£] given by eqn. (|3 .39[) as 

Since the noise probability distribution function is Gaussian the only necessary correlation functions for the 
noise are given by 

«Mt))> = 0, (<£„,£(*) W (*')>) =>C a p(k;t-t')6 3 (k + k') (3.45) 

and the higher order correlation functions are obtained from Wick's theorem as befits a Gaussian distribution 
function. Because the noise kernel JC a p(k; t — t') ^ 8(t — t') the noise is colored. 

• The average over the initial conditions with the Wigner distribution function W(^E' , II ) which we denote as 

jD^°jDn° w(* ;n°)^[*°,n ] 



" 4[ * ' n ° ] " JD^JDW w(*° ; n°) (3 - 46) 

Therefore, the average in the time evolved reduced density matrix implies two distinct averages: an average over 
the initial conditions of the system fields and and average over the noise distribution function. The total average is 
defined by 



(...) = «...)) . (3.47) 

Equal time expectation values and correlation functions are simply expressed in terms of the center of mass Wigner 
variable iff as can be seen as follows: the expectation values of the field 

(<t> + (x,t))=Tr<t>(x,t)p(t) ■ (4>-(x,t)) = Trp(t)<t>(x,t) (3.48) 
hence the total average (|3.47p is given by 

(<j>(x,t)) = (V(x,t)). (3.49) 
Similarly the equal time correlation functions obey 

{<j> + {x,t)<j> + {x",t)) = (<t> + (2,t)4r{Z,t)) = (<t>-(x-,t)<p+(x',t)) = {4>-@,t)4>-(#,t)) = (*(2,t)9(3?,t)). (3.50) 

Therefore the center of mass variables iff contain all the information necessary to obtain expectation values and equal 
time correlation functions. 
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A. One body density matrix and equilibration: 

We study equilibration by focusing on the one-body density matrix 

PaP {k-t) = Trp(o)<£ a (M)M-M) = H7>(t)0a(*,o)M-M) 



where 



p(t) = e- lHt p(0)e lHt 



(3.51) 



(3.52) 



is the time evolved density matrix. The time evolution of the one-body density matrix obeys the Liouville-type 
equation 



■£Pap(t) = -iTr[H, p{t)} <j> a (k, 0)M-k, 0) . 
If the system reaches equilibrium with the bath at long times, then it is expected that 



(3.53) 



[H, P (t)] 







(3.54) 



Therefore the asymptotically long time limit of the one-body density matrix yields information on whether the density 
matrix is diagonal in the flavor or any other basis. Hence we seek to obtain 



Pa p{k- oo) = TVp(oc)<£ a (fe, 0)^/3 (-M) = a% {oo)^ _ % {co)) 



(3.55) 



and to establish the basis in which it is nearly diagonal. The second equality in eqn. (|3.55[) follows from eq. (|3.50[) . 
and the average is defined by eq. (|3.47p . To establish a guide post, consider the one-body density matrix for the free 
"neutrino fields" in thermal equilibrium, for which the equilibrium density matrix is 



n - p -/5HoM 



(3.56) 



where Hq [tp] is the free "neutrino" Hamiltonian. This density matrix is diagonal in the basis of mass eigenstates and 
so is the one-body density matrix which in the mass basis is given by 



Pij(k) = 



2uj 1 (k) 



coth 




2 





2^(fcy C0th 



2 



(3.57) 



therefore in the flavor basis the one-body density matrix is given by 





p a0 {k) = U{6) 



2^7W coth 




^i(fc) 

2 



2w 2 (k) 



coth 



2 



(3.58) 



This simple example provides a guide to interpret the approach to equilibrium. Including interactions there is a 
competition between the mass and flavor basis. The interaction is diagonal in the flavor basis, while the unperturbed 
Hamiltonian is diagonal in the mass basis, this of course is the main physical reason behind neutrino oscillations. In 
the presence of interactions, the correct form of the equilibrium one-body density matrix can only be obtained from 
the asymptotic long time limit of the time-evolved density matrix. 



B. Generalized fluctuation-dissipation relation: 

From the expressions (|3.35[) and (|3.36[) in terms of the Wightmann functions (|3.20I3.21[) which are averages in the 
equilibrium density matrix of the bath fields (%, W), we obtain a dispersive representation for the kernels JC a /3(k; t — 
t');T,^ /3 (k;t — t'). This is achieved by writing the expectation value in terms of energy eigenstates of the bath, 
introducing the identity in this basis, and using the time evolution of the Heisenberg field operators to obtain 

/oo />oo 
dej a> (k,u) e^*"*') ; G 2 Q^{k;t-t') = / du; a< (k,oj) e 1 ^'^ (3.59) 
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with the spectral functions 



G 2 



< (k,u) = —J2^ En (MO a ^rn){m\O^(0)\n)5(uj-(E n -E m )) 

m.n 

= e- pEm {n\G ^ k mm){m\O a ^)\n) 6(u - (E m - E n )) 



(3.60) 
(3.61) 



where Zf, = Tr e~@ Hx is the equilibrium partition function of the "bath" and in the above expressions the averages 
are solely with respect to the bath variables. Upon relabelling m <-> n in the sum in the definition (|3 . 6 1 [) and using 
the fact that these correlation functions are parity and rotational invariant [57| and diagonal in the flavor basis we find 
the Kubo-Martin-Schwinger (KMS) relation \El\ 

<r^(k,w) = cr> ap {k-u) = e^a> (k,u) . (3.62) 

Using the spectral representation of the 8(i — t') we find the following representation for the retarded self-energy 



with 



^(fc;fco) 



( 

2tt 



*° £* (k; ho) 



[°'«/3( fc ; w )- cr a/3( fc ; w )] 



u> — ko + ie 

Using the condition (|3.62p the above spectral representation can be written in a more useful manner as 



£&(*;*o) 



ImE^a/) 



dio- 

^ uj-ko + ie 



where the imaginary part of the self-energy is given by 

ImE* = 7ra>^k;cj)[e^-l\ 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



and is positive for u) > 0. Equation (|3.62p entails that the imaginary part of the retarded self-energy is an odd function 
of frequency, namely 



ImE^(fe;a;) = -ImE^(fe;-a;) . 
The relation (|3. 66[) leads to the following results which will be useful later 



<li(k-u) = -lmY,%{k;u>)n{u) ; <r<,(fc;w) = -ImE^(fc;o;) [1 + n(w)] 



(3.67) 



(3.68) 



where n{uj) — \e^ u — 1] 1 is the Bose-Einstein distribution function. Similarly from the definitions (|3.35[) and (|3.59|) 
and the condition (|3.62p we find 



K af} (k;t-t') = /" ^e^C-'U a0 (k;k o ) 
J-oo 2tt 

iC af3 (k;k ) = 7ra>p{k;k ) [e 0k ° + l] 
whereupon using the condition (13.62p leads to the generalized fluctuation-dissipation relation 

\0ko~ 



/C Q(3 (fc; k ) = ImE Q/3 (fc; k ) coth 



Thus we see that ImE^^fc; k Q ) ; K. a0 (k; k ) are odd and even functions of frequency respectively. 
For the analysis below we also need the following representation (see eqn. (|3.36p ) 



(3.69) 
(3.70) 

(3.71) 
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E a/ ,(*;t-t')=-t / e^-V 

J — oo 

whose Laplace transform is given by 



do; 



e ^(*-*')l mS fl Jk;u)dw (3.72) 



1 /■<» ImE^u;) 



E a/3 (fc;s)= / dte- st Z al3 (k;t) = 

This spectral representation, combined with (|3 . 65[) lead to the relation 

fco) = E a/3 (fc; s = ifc + e) 



(3.73) 



(3.74) 



The self energy and noise correlation kernels E,/C are diagonal in the flavor basis because the interaction is diagonal 
in this basis. Namely, in the flavor basis 



E(fc,w) 



S ee (fc,o;) 
S MM (fc,w) 



/C = [l + 2n(w)]ImE(A;,a;) 



K, ee (k,uj) 
/C MM (fc,a;) 



In the mass basis these kernels are given by 



E = -(E ee + E w ] 11 



v ~ \ / cos 29 sin 29 
Lee "M sin 26» -cos26> 



and 



(3.75) 



(3.76) 



= - ( ^ee + ^ W I 1 



cos 26» sin 26» 
w ; I sin 20 -cos 26* 



(3.77) 



IV. DYNAMICS: SOLVING THE LANGEVIN EQUATION 

The solution of the equation of motion (|3.42p can be found by Laplace transform. Define the Laplace transforms 



(4.1) 



along with the Laplace transform of the self-energy given by eqn. (|3.73p . The Langevin equation in Laplace variable 
becomes the following algebraic matrix equation 



(s 2 + k 2 )8 aP + M 2 a0 + Y af} + E a/3 (fc; s) 



(4.2) 



where we have used the initial conditions (|3 .41 [) . The solution in real time can be written in a more compact manner 
as follows. Introduce the matrix function 



G(k;s) 



(s z + r)l + M + ¥ + E(fc; s) 



and its anti-Laplace transform 



f ds ~ 



(4.3) 



(4.4) 
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where C refers to the Bromwich contour, parallel to the imaginary axis in the complex s plane to the right of all the 
singularities of G(k;s). This function obeys the initial conditions 



G Q/3 (fc;0)=0 ; G Q/3 (fc;0) = l. 
In terms of this auxiliary function the solution of the Langevin equation (|3.42[) in real time is given by 



y a i (t; *°; n°; = G a p (*! *) + G a p (fc; t) n° % + / G aP (k; t') ^(i - t')dt' , 

« o 

where the dot stands for derivative with respect to time. In the flavor basis we find 

SM 2 ( cos26>- A(fc;s) -sin 29 



G f {k;s) = S{k;s) 



whereas in the mass basis we find 



+ ZJ 2 (k) + S(fc;s) 1 + 



sin 29 



cos 26* + A(k;s) 



G m (k; s) = S(k; s) 



(,s 2 +cJ 2 (fc) + S(fc;.s))l + 



SM 2 I 1 - A(fc;s) cos 20 A(fc;s)sin20 
2 I A(fc;s)sin26» -1 + A(fc; s) cos 29 



where M and 5M 2 were defined by eqn. (|2.7[) and to simplify notation we defined 



ZJ(k) = \/k 2 + M 
S(fc;s) = -(£ee(fc;s)+V r ee +£ w (A;;s)+V r w ) 



A(k;s) 



M 2 - M 2 



and 



S{k;s) 



2 r ii jf 2 2 

(s 2 +LJ 2 (k) + £(*;*)) - (— — ) [(cos 26* - A(fc;s)) 2 + (sin26») 2 



(4.5) 



(4.6) 



(4.7) 



(4.E 



(4.9) 
(4.10) 

(4.11) 



(4.12) 



In what follows we define the analytic continuation of the quantities defined above with the same nomenclature to 
avoid introducing further notation, namely 



uj) = s = ilo + e) ; A(k; u>) = A(fc; s = ilo + e) . 
Their real and imaginary parts are given by 

T, R (k;uj) = - [E_R jee (fc, ui) + S (k,u) + V ee + V^] 



A a (*;w) 
A/(fc;w) 



1 



<5M 2 
1 



ttt [Sfl,ee(fc, w) - S RiAlAl (A;,a;) + V ee - F A 



<5M 



(4.13) 

(4.14) 
(4.15) 
(4.16) 
(4.17) 



We remark that while the matter potential V is of of order G, X is of order G 2 . Therefore, in perturbation theory 

» Zi(k;u) ; A R (k; u) » A/(A; w) . (4.18) 

This inequality also holds in the standard model, where the matter potential is of order Gf[4^| while the absorptive 
part that determines the relaxation rates is of order G 2 F . This perturbative inequality will be used repeatedly in the 
analysis that follows, and we emphasize that it holds in the correct description of neutrinos propagating in a medium. 
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V. SINGLE PARTICLES AND QUASIPARTICLES 

Exact single particle states are determined by the position of the isolated poles of the Green's function in the 
complex s— plane. Before we study the interacting case, it proves illuminating to first study the free, non-interacting 
case. 



A. Free case: G — 

To begin the analysis, an example helps to clarify this formulation: consider the non-interacting case G = in 
which £ = 0; A = 0. In this case G/ jm (fc; s) have simple poles at s = ±iu>i(k) and ±iu>2{k) where 



Wi(k) = ^k 2 + Mf ; i= 1,2. 
Computing the residues at these simple poles we find in the flavor basis 

Gf (jb; t) = BinMfc)*) E(1 ) {e) + 4f) R ( 2 ) 

where we have introduced the matrices 

'l o' 



R (1) (6>) 



cos a 
cos 9 sin ( 



cos V sin f 
sin 2 9 



U(0) I o J U ~ X W 



(5.1) 



(5.2) 



(5.3) 



In the mass basis 



with the relation 



R( 2 \9) = 



Sin 9 cosflsinfl \ x 

cos 9 sin 9 cos z 9 / w 1 w 



G m (k: t) — 



sin(uJi (k)t) 
^l(fc) 





s'm(uJ2(k)t) 



G f {k;t) = U{9) G m (k;t) U~\ff) 



(5.4) 



(5.5) 



(5.6) 



and U(9) is given by 12.41 Consider for simplicity an initial condition with ^° ^ 0; II = in both cases, flavor and 
mass. The expectation value of the flavor fields $ a in the reduced density matrix (|3.47p is given by 



and that for the fields in the mass basis is 



cos 



(cJi(k)t) RW(6) + cos(uj 2 (k)t) RW(9) 



e,k 



*!(M)\\ _ ( *ifc C0S K( fc )*) \ 

V 3 (k;t))/ ~ ^ tt^cosM*:)*) ) 



(5.7) 



(5.8) 



These are precisely the solutions of the classical equations of motion in terms of flavor and mass eigenstates, namely 

<j> e (k;t) = cos 9 (pi (k; t) + sin 9 <p2 (k; t) 

</> M (fc;i) = cos9 tf2(k;t) — sin9 (f2{k;t) (5.9) 

where 



<Pi{k; t) = <pi(k; 0) coswi(fc)i ; (f2{k; t) = <p 2 (k; 0) cos U2(k)t 



(5.10) 
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for vanishing initial canonical momentum and the initial values are given in terms of flavor fields by 



(pi(k;0) — cos0 e (fc;O) — sin6'0 / i(A;;O) 
<p 2 (k;0) = cos0 M (fc; 0) + sin0^ e (fc;O) 



(5.11) 



inserting (|5. 10|) with the initial conditions (|5.11|) one recognizes that the solution (15. 7p is the expectation value of the 
classical equation of motion with initial conditions on the flavor fields and vanishing initial canonical momentum. 

It is convenient to separate the positive (particles) and negative (antiparticles) frequency components by considering 
an initial condition with II 7^ 0, in such a way that the time dependence is determined by phases corresponding 

to the propagation of particles (or antiparticles). Without mixing (9 = 0) this is achieved by choosing the following 
initial conditions 

n° r = ^irtJk) *° r (5.12) 

a.k ' x ' a.k v ' 

for particles (— ) and antiparticles (+) respectively as in eq. (|2.11| . This choice of initial conditions leads to the result 



cos(wi (k)t) =F i—^—TTT sa\(uoi{k)t) 



cos(u> 2 (fc)i) T i ^ ^ sm(Lo 2 {k)t) 
uj 2 {k) 



(5-13) 

It is clear from (|5 . 1 3|) that no single choice of frequencies £lp(k) can lead uniquely to time evolution in terms of 
single particle/ antiparticle phases e T^i,2(fc)* p This is a consequence of the ambiguity in the definition of flavor states 
as discussed in detail in section (|II Aj) . However, for the cases in which \M 2 — M|| <C (k 2 + M ) , relevant for 

relativistic mixed neutrinos, and for K°K and B°B mixing, the positive and negative frequency components can 
be approximately projected out as follows. Define 



U(k) = \lk 2 + M 

in the nearly degenerate or relativistic regime when \8M 2 \/uj 2 (k) < 1 



(5.14) 



u>i(k) = ui{k) 



SM 2 I SM 2 



4ZJ 2 (k) \uJ 2 (k) 



; uj 2 (k) = uj(k) 



SM 2 ( 5M 2 



4ZJ 2 (k) \uJ 2 (k) 



(5.15) 



Taking Clp(k) = u>(k) and choosing for example the negative sign (positive frequency component) in (|5.13|) we find 



SM 2 



Consider the following initial condition 



neglecting the corrections in (|5.16p we find 



= e,fc 

k \ 



(5.16) 



(5.17) 



SM 2 



cJ 2 (fc) 



(5.18) 



which is identified with the usual result for the oscillation transition probability ^> e — > Vf^ upon neglecting the 
corrections. 
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B. Interacting theory, G^O 

For G 7^ 0, the self-energy as a function of frequency and momentum is in general complex, the imaginary part 
arises from multiparticle thresholds. When the imaginary part of the self-energy does not vanish at the value of the 
frequency corresponding to the dispersion relation of the free particle states, these particles can decay and no longer 
appear as asymptotic states. The poles in the Green's function move off the physical sheet into a higher Riemann 
sheet, the particles now become resonances. 

Single particle states correspond to true poles of the propagator (Green's function) in the physical sheet, which are 
necessarily away from the multiparticle thresholds. This case is depicted in fig. 



FIG. 2: Bromwich contour in s-plane, the shaded region denotes the cut discontinuity from multiparticle thresholds across the 
imaginary axis, the filled circles represent the single particle poles. 



Let us consider the Green's function in the flavor basis eqn. (|4.7[) . The single particle poles are determined by the 
poles of S(k; s) on the imaginary axis away from the multiparticle cuts. These are determined by the roots of the 
following equations 



n?(fc) = Zj 2 (fc) + £fl(fc;fti(fc))- 

nl(fc) = Lj 2 (k) + z R (k;n 2 (k)) + 

along with the conditions 



SM 2 
2 

SM 2 



cos 20- A fl (fc;fti(fc))) : 
^ cos 2(9- A R (k;Cl 2 {k)))' 



sin 26 



sin 26 



£i(fc;Qi, 2 (fc)) = ; A/(fc;Jli, 2 (fc))=0 



(5.19) 
(5.20) 

(5.21) 



where the subscripts -R, / refer to the real and imaginary parts respectively. Evaluating the residues at the single 
particle poles and using that the real and imaginary parts of the self-energies are even and odd functions of frequency 
respectively, we find 



(5.22) 



where Gf yCU t(t) is the contribution from the multiparticle cut, the matrices are given by (|5.3I5.4|) and 9^ 2 (k) 

are the mixing angles in the medium 



cos2^(fc) 



cos 20- A R (Qi(k)) 



'cos2d- A R (k;ni(k))) 2 + (sin20) 2 



- ; sin 20^*;) 



sin 20 



(cos 20- A R (k ; n t (k))) 



2 + (sm28) 2 



for i = 1,2. The wave function rcnormalization constants are given by 



(5.23) 
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cos2C(fc)Afl(fc;w) 



(5.24) 



=n*(fc) 



where the prime stands for derivative with respect to w. At asymptotically long time the contribution from the cut 
Gf,cut(t) ~ t~ a where a is determined by the behavior of the self-energy at threshold(55l. |5(| . 

In perturbation theory and in the limit W{k) 2 \SM 2 | the dispersion relations (|5. 19I5.20[) can be solved by writing 



±n l (k) = ±(u(k) + Sun(k)), 



(5.25) 



we find 



where we defined 



2u>{k) 4w(k) 



(5.26) 



g(k;uj) 



and the shorthand 



(cos 20 - A R (k; w)f + (sin 2Bf 



g(k) = g(k;uj = u>(k)) 



(5.27) 



(5.28) 



To leading order in the perturbative expansion and in 5M 2 /u> 2 (k) we find 6$(k) = ( £\k) = 6 m (k). Gathering these 
results, we find the dispersion relations and mixing angles in the medium to be given by the following relations, 



Ui(k) = u>(k) H — — g(fc) 

2 57(A) 457(fc) 

r, ,n S R (fc;5J(fc)) <5M 2 .,, 



and 

cos 29 m {k) 



cos 26* - A R (k;uJ(k)) 



( cos26 ~ A R {k;ZJ(k))) 



2 + (sin 2(9) 2 



25J(fc) 4ZZJ(fc) 



— ; sin20 m (£;) 



(5.29) 
(5.30) 



sin 20 



(cos 20- Ajj(fc;aT(fc)))' 



sin 



20)' 



(5-31) 

These dispersion relations and mixing angles have exactly the same form as those obtained in the field theoretical 
studies of neutrino mixing in a medium (49l l54j . 



C. Quasiparticles and relaxation. 



Even a particle that is stable in the vacuum acquires a width in the medium through several processes, such as 
collisional broadening or Landau damping In this case there are no isolated poles in the Green's function in the 
physical sheet, the poles move off the imaginary axis in the complex s— plane on to a second or higher Riemann sheet. 
The Green's function now features branch cut discontinuities across the imaginary axis perhaps with isolated regions 
of analyticity. The inverse Laplace transform is now carried out by wrapping around the imaginary axis as shown in 
fig. (O, and the real time Green's function is given by 

G a p(k;t) = - duje lLJt lmG a0 (k;s = ico + e) (5.32) 

T J-oc 

Under the validity of perturbation theory, when the inequality (|4.18p is fulfilled we consistently keep terms up to 
0(G 2 ) and find the imaginary part to be given by the following expression 



(5.33) 
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where we have introduced 



D±(k;w) = -lj 2 +ujI + Y lR {k;uj) T^SM 2 g(k;Lu) 

y±(k;ui) = - (1 ± cos20 m (w, k)) S/, ee (fc; + - (1 T cos20 m (w, k)) Ej jWt (fc;o;) 

are the real and imaginary parts of the self energy respectively, with 

1 



A R (k;uj) = \E R , ee (k;io) + V ee - T, Rtfifi (k;w) - V w ] , 



and 



A(k; u) = [~lu 2 + Zuj 2 + Y, R (k; u)l] + 



SM 2 ( cas26- A R (k;u) 



sin 29 



-sin 26 -cos 26* + A R (k;u;) 



(5.34) 
(5.35) 

(5.36) 

(5.37) 



The denominator in (|5.33|) features complex zeroes for 

D ± {k-u J )+ 1± {k-u J )=Q. 



(5.38) 



(5.39) 



Near these zeroes lmG(k; s = iu> + e) has the typical Breit-Wigner form for resonances. The dynamical evolution 
at long times is dominated by the complex poles in the upper half lu— plane associated with these resonances. In 
perturbation theory the complex poles of ImG(fc; s = iu; + e) occur at 



and 



lu = ±Qi(k) + i 



lu = ±Q 2 (k) + i 



ri(fc) 



r 2 (fc) 



where £lx t 2{k) are given by (I5.19|5.20| and 

r^fc) _ 7+ (fe s n 1 (A)) _ 7+(*.w(*)) . r 2 (fc) _ 7 _(fc,r! 2 (fc)) _ 7 _(fc,cj(fc)) 



2 2fii(fc) 257(A) ' 2 

These relaxation rates can be written in an illuminating manner 



2n 2 (k) 



2ou(k) 



(5.40) 



(5.41) 



(5.42) 



ri(fc) = r ee (fc)cos 2 Mfc) + r w (fc)sin 2 Mfc) 

T 2 (fc) = T^(k) cos 2 6 m (k) + T ee (k) sin 2 m (k) 



where 



-Tact \k) 



Ej, aa (fc;u;(fc)) 
Zu(k) 



(5.43) 
(5.44) 



(5.45) 



are the relaxation rates of the flavor fields in absence of mixing. These relaxation rates are similar to those proposed 
within the context of flavor conversions in supernovae[58|, or active-sterile oscillations [28l I50ll59j. 

We carry out the frequency integral in (|5.32[) by approximating the integrand as a sum of two Breit-Wigner 
Lorentzians near lu = ±£2^2 (&) with the following result in the flavor basis, 



G f (k;t)= Z^e-^ 



\li[K) 1 \l\[k) 

sin(0 2 (fc)t) E ( 2)( ^ (fe)) + Wl coBiMQt) R (3) (0m(fc)) 



n 2 (fc) 



n 2 (k) 



(5.46) 
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FIG. 3: Bromwich contour in s-plane, the shaded region denotes the cut discontinuity from multiparticle thresholds across the 
imaginary axis. 



where again we have used the approximation \SM \ -C lu (k) and introduced 



^v. = si„2, ( ±Z ™l 2 l)=^°m (11) u-w 



and 



7(A) = 



cos 26» - sin 20 / 

E/ iee (fc;uJ(fc) - £r lWt (fe; 57(A)) 



SM 2 g(k) 



(5.47) 



(5.48) 



Under the assumption that S_r jQ , q ^> £/, Q .a it follows that 7(fe) < 1. As in the previous section, we can approxi- 
mately project the positive frequency component by choosing the initial condition (|5.12j) with il a = 57(fc), leading to 
the result 



<*.,«(*)> 



f'e^'*-^ (R (1 HM*0)+i^R (3) (M*))) + 



^)ku(^) 

u 2 (fc) / J v 



where 



W(fc) = uJ(fc) 

r(A) i 



2 

Aw(fc) 

Ar(fc) 



4ZJ(fc) 
2 £(fc) 
4w(fc) 
cos 20 m 
4ZJ(fc) 



S fl (fc;57(fc)) 
4u7(fc) 

S/ iee (fc,ZU(fc) + Ej )W ,(A;,5J(fc) 



,(k,oj(k) - T, Itflfl (k,u!(k) 



With the initial condition (|5.17p we now find the long time evolution of the transition probability iff e — > 



2 sin 2 29 m (k) 



-rxWt + e -r 3 (fc)t _ 2 e -i(r l(fc )+r 2 ( fc ))t cog ( 2Aw (k)t) 



e.A; 



(5.50) 
(5.51) 
(5.52) 
(5.53) 

(5.54) 



where we have neglected perturbatively small terms by setting ~ 1 ; j(k) ~ in prefactors. The solution (|5.49D 
can be written in the following illuminating form 



U- X {9 m {k)) & % (5.55) 
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where 



n(o fiA\ - ( cos9 m (k)(l + ij(k)) ain.0 m (k)(l - vy{k)) 
M{u m (K)j — i - S ine m (k) cos 6 m {k) 



(5.56) 



7 J- 1 (ft (hW - - ( cose m{k) -sm6 m (k)(l-ij(k))\ , . 

lroljj " + i C o S 6 m (k)j(k))\-sme m (k) cos0 m (*)(l+*7(*)) J 1 j 

Obviously the matrix is not unitary. 

D. Long time dynamics: Weisskopf-Wigner Hamiltonian 

A phenomenological description of the dynamics of mixing and decay for neutral flavored mesons, for example 
KqKo ; -Bo-Bo relies on the Weisskopf-Wigner (WW) approximation [60|. In this approximation the time evolution 
of states is determined by a non-hermitian Hamiltonian that includes in a phenomenological manner the exponential 
relaxation associated with the decay of the neutral mesons. This approach has received revived attention recently 
with the possibility of observation of quantum mechanical coherence effects, in particular CP-violation in current 
experiments with neutral mesons [6 lj. In ref.[62j a field theoretical analysis of the (WW) approximation has been 
provided for the K Ko system. 

The form of the solution (|5.55| suggests that a (WW) approximate description of the asymptotic dynamics in 
terms of a non-hermitian Hamiltonian is available. To achieve this formulation we separate explicitly the fast time 
dependence via the phase e TM fe )* for the positive and negative frequency components, writing 

*g(t) = e ~^( fc )* $+( t ) + e l ~( fe )* *-( t ) (5.58) 



where \&g (t) the amplitudes of the flavor vectors that evolve slowly in time. The solution for the positive frequency 



k 

component (|5.49[) follows from the time evolution of the slow component determined by 

d 



dt 

where the effective non-hermitian Hamiltonian TL is given by 

SAP /-cos26> sin2(9\ , aJ(fc)) + V . £ 7 (fc;57(fc)) 



XM- (-COS26 sin 20 ^ , ^r{K} ) + v ■ ^A fc i ^(k)) , . 

" AU(k) \ sin 20 cos2^ + 2 U(fc) * 2tJ(fc) 1 j 



with 



+ * - + K„ Wfc ,» W) + Vw ) (5.61, 

The (WW) Hamiltonian 7^ can be written as 

_ f S K (fc;^(fc)) J(fc)^ o~Af 2 g(k) ( -(cos20 m (k) + ij(k)) sin2^ m (fc) \ 

4U(k) 1 2 J + Au{k) \ sm29 m (k) (cos 20 m (fc) + ij(k)) J 

where we have used the definitions given in equations (|5. 2715. 2815. 3115. 4815. 51|) . It is straightforward to confirm that 
the Wigner-Weisskopf Hamiltonian can be written as 

H w =U{B m {k)) I A - Q (fc) A+ ° (fc) J U^iOmik)) (5.64) 
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where lA(9 m {k)) is given in (|5.56|) and using the definitions given in eqn. (j5.50H5.53f the complex eigenvalues are 

A T (fc) = W(k) - U{k) - T yAu(k) + i^j^-J (5.65) 



The solution of the effective equation for the slow amplitudes (I5.59|) coincides with the long time dynamics given 
by (|5.49p when the wave function renormalization constants are approximated as Z^S 1 ~ 1. Therefore the (WW) 
description of the time evolution based on the non-hermitian Hamiltonian TL W (|5.60p effectively describes the evolution 
of flavor multiplets under the following approximations: 

• Only the long-time dynamics can be extracted from the Weisskopf-Wigner Hamiltonian. 

• The validity of the perturbative expansion, and of the condition SM 2 <C ZU(fc) 2 . 

• Wavefunction renormalization corrections are neglected Z^ ~ 1 and only leading order corrections of order 
7(fc) are included. 

While the Weisskopf-Wigner effective description describes the relaxation of the flavor fields, it misses the stochastic 
noise from the bath, therefore, it does not reliably describe the approach to equilibrium. 

VI. EQUILIBRATION: EFFECTIVE HAMILTONIAN IN THE MEDIUM 



As discussed in section llll Al we study equilibration by focusing on the asymptotic long time behavior of the one-body 
density matrix or equal time correlation function, namely 



,!™<*«,s(*)*0,-i;(*)>- 



(6.1) 



In particular we seek to understand which basis diagonalizes the equilibrium density matrix. 

Consider general initial conditions ^° ^ and II ^ 0, in which case the flavor field *Sf z(t) is given by Eq. 

(|4.6j) with Gf(k;t) given by eqn. (I5.46| . For t IY2, the first two contributions to (|4.6|) which depend on the 

initial conditions fall-off exponentially as e 2 and only the last term, the convolution with the noise, survives at 
asymptotically long time, indicating that the equilibrium state is insensitive to the initial conditions as it must be. 

To leading order in the perturbative expansion in G, and in the limit SM 2 /uJ 2 (k) <C 1, we can approximate 
9rn (k) ss 9m\k) = 9 m (k), where the effective mixing angle in the medium 9 rn (k) is determined by the relations 
(I5.31[) . Similarly we can approximate the wave function renormalization constants as Z^ x \k) Z^ 2 \k) = Z(k) with 



Z(k) = 



1 - -^(z R (k;Lu) - i-lY^f- cos2e m (k)A R (k; u) 



1 -1 



-1 LJ-Lu(k) 

where the prime stands for derivative with respect to lu. Thus, Gf{k;f) and G m (k;t) are related by 

G f (k; t) » Z(k) U(9 m ) G m (k; t) U-\9 m ) 

where G m (k;t) is given by 



G m (k;t) 





sin(n 2 (fc)t) -EaMt 
o 2 (fe) e 



4^ sin20 m (fc) 



r 2 CO t COslQj (k)t) rjXjQ f COsf^i (k)t) 



fi 2 (fc) 









1) 



It is useful to define the quantities h m (t,uj) and £^ as follows 



(6.2) 



(6.3) 



(6.4) 



h m (t,uj) = / e- iut ' G m (k;t)dt' 
Jo 



(6.5) 
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and 

r+oo 



fait -t')= e^C^^gMdw, (6.6) 

J — oo 

with the noise average in the flavor basis given by 

(( f Pl J?M = ^(A; + w') = ImS« (A; u) coth f^) % + u') . (6.7) 

We find convenient to introduce 

K m (k;uj) = U- 1 (e m )jC(k;u)U(6 m ) . (6.8) 

The approach to equilibrium for t » T^J, can be established from the unequal time two-point correlation function, 
given by 

Jjm^ a % (t)^ _ u (t'))=Z\k)U(9 m )S^J + du e^-*') h m (oo,co) ]C m (k; u) h m (oo, -oojtr 1 ^) (6.9) 

where we have taken the upper limit t — ► oo in (|6.6[) . The fact that the correlation function becomes a function of 
the time difference, namely time translational invariant, indicates that the density matrix commutes with the total 
Hamiltonian in the long time limit. The one-body density matrix is obtained from (|6.9|) in the coincidence limit t = t' . 
Performing the integration over u>, we obtain after a lengthy but straightforward calculation 

& (*a,5(*)*A-5(*)> - V<M ( tS) tS) ) Erl(ffw)5 (6 ' 10) 



wherein 



A ll( fc) = coth f «2 ) . A22(fc) = _L_ coth f ] , (6.11) 

v 7 20i(fc) V 2 / 20 2 (fc) V 2 j 



and to the leading order of SM 2 /u>(k) -C 1, we find A21 = Ai2(f2i — > ^2) where 

a ™ 1 „f ^( k ) \ ■ on WhM 2n 1 (fc)(r 1 (fc) + r 2 (fc)) 

Al2(fc) = 27Mfc) coth (-2- J sin2 ^ (fc) i + (,(fc)) 2 ; ^ = 5M^(fe) ■ (6 - 12) 

Since 7(A) -C 1, it is obvious that Ai2(A) and A2i(fc) are perturbatively small compared with An(k) and A22(A), 
in either case 77(A) S> 1 or n(k) <C 1. The asymptotic one-body density matrix (|3.55p then becomes 



z 



■ coth 



^Qi(fc) ^ 



p a ,p(k;oc) = U(9 m ) I 2 ^W^V 2 y !/ fl0 , M N I t/" 1 ^) ; e<0(G 2 ) (6.13) 



z 

2n 2 (k) 



coth 



where we neglected corrections of 0(G 2 ) in the diagonal matrix elements. 

Neglecting the perturbative off-diagonal corrections, the one-body density matrix commutes with the effective 
Hamiltonian in the medium which in the flavor basis is given by 

H ff (k)-u(k)( 1 0ViW-cos2fl wx20\ ._}_( Z R Mk-Mk)) + V ee \ 

this effective in-medium Hamiltonian can be written in a more illuminating form 

H eff (k) = U(9 m ) ( j U-\e m ) (6.15) 
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where f2i,2(fc) are the correct propagation frequencies in the medium given by eqn. (|5.29I5.30|) . 

This effective Hamiltonian includes the radiative corrections in the medium via the flavor diagonal self-energies 
(forward scattering) and apart from the term proportional to the identity is identified with the real part of the 
Weisskopf-Wigner Hamiltonian Ti w given by eqn. (|5.60[) . This form highlights that the off-diagonal elements of the 
one-body density matrix in the basis of eigenstates of the effective Hamiltonian in the medium are perturbatively 
small. The unitary transformation U(0 m ) relates the flavor fields to the fields in the basis of the effective Hamiltonian 
in the medium. 

Comparing this result to the free field case in thermal equilibrium, where the one body density matrix in the flavor 
basis is given by eqn. I|3.58[) . it becomes clear that in the long time limit equilibration is achieved and the one-body 
density matrix is nearly diagonal in the basis of the eigenstates of the effective Hamiltonian in the medium (|6.14p 
with the diagonal elements determined by the distribution function of these eigenstates. 

This means that within the realm of validity of perturbation theory, the equilibrium correlation function is nearly 
diag onal in the basis of the effective Hamiltonian in the medium. This result confirms the arguments advanced in 
4f]. Since the effective action is quadratic in the "neutrino fields" higher correlation functions are obtained as Wick 
contractions of the two point correlators, hence the fact that the two point correlation function and consequently 
the one-body density matrix are diagonal in the basis of the eigenstates of the effective Hamiltonian in the medium 
guarantee that all higher correlation functions are also diagonal in this basis. 

A. On "sterile neutrinos" 

The results obtained in the previous sections apply to the case of two "flavored neutrinos" both in interaction with 
the bath. However, these results can be simply extrapolated to the case of one "active" and one "sterile" neutrino 
that mix via a mass matrix that is off-diagonal in the flavor basis. By definition a "sterile" neutrino does not interact 
with hadrons, quarks or charged leptons, therefore for this species there are no radiative corrections. Consider for 
example that the "muon neutrino" represented by M does not couple to the bath, but it does couple to the "electron 
neutrino" solely through the mixing in the mass matrix. Since the interaction is diagonal in the flavor basis, the 
decoupling of this "sterile neutrino" can be accounted for simply by imposing the following "sterility conditions" for 
the matter potential Y and the self energies 

= ; S fl , w = ; E JiW , = . (6.16) 

All of the results obtained above for the dispersion relations and relaxation rates apply to this case by simply imposing 
these "sterility conditions" . In particular it follows that 

rx(fc) = r ee (fc) cos 2 8 m (k) ; r 2 (fc) = r ee sin 2 e m (k) (6.17) 

where T ee (k) is the relaxation rate of the active neutrino in absence of mixing. This result highlights that in the 
limit 8 — ► the in-medium eigenstate labeled "2" is seen to correspond to the sterile state, because in the absence of 
mixing this state does not acquire a width. However, for non-vanishing vacuum mixing angle, the "sterile neutrino" 
nonetheless equilibrates with the bath as a consequence of the "active-sterile" mixing, which effectively induces a 
coupling between the "sterile" and the bath[27], HH, [H,[5!|. The result for ^(fc), namely the relaxation rate of the 
"sterile" neutrino is of the same form as that proposed in refs. [13, HH, [5(J, [H, [59[ . The result for the "sterile" rate ^(fc) 
compares to those in these references in the limit in which perturbation theory is valid, namely S ee (fc) /SM 2 ~g(k) <C 1 
since the denominator in this ratio is proportional to the oscillation frequency in the medium. 

VII. SUMMARY OF RESULTS AND CONCLUSIONS 

In this article we studied the non-equilibrium dynamics of mixing, oscillations and equilibration in a model field 
theory that bears all of the relevant features of the standard model of neutrinos augmented by a mass matrix off 
diagonal in the flavor basis. To avoid the complications associated with the spinor nature of the neutrino fields, we 
studied an interacting model of flavored neutral mesons. Two species of flavored neutral mesons play the role of two 
flavors of neutrinos, these are coupled to other mesons which play the role of hadrons or quarks and charged leptons, 
via flavor diagonal interactions that model charged currents in the standard model. These latter meson fields are 
taken to describe a bath in thermal equilibrium, and the meson-neutrino fields are taken to be the "system" . We 
obtain a reduced density matrix and the non-equilibrium effective action for the "neutrinos" by integrating out the 
bath degrees of freedom up to second order in the coupling in the full time-evolved density matrix. 
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The non-equilibrium effective action yields all the information on the particle and quasiparticle modes in the 
medium, and the approach to equilibrium. 
Summary of results: 

• We obtain the dispersion relations, mixing angles and relaxation rates of the two quasiparticle modes in the 
medium. The dispersion relations and mixing angles are of the same form as those obtained for neutrinos in a 
medium[ii,[53. 

• The relaxation rates are found to be 

ri(fc) = T ee (k)cos 2 e m (k) + r^(k)sm 2 e m (k) (7.1) 
r 2 (fc) = r^{k)cos 2 e m (k) + r ee {k)sm 2 e m (k) (7.2) 

where 

LU(k) 

are the relaxation rates of the flavor fields in absence of mixing and ~Ei, a ,a are the imaginary parts of the 
"neutrino" self ener gy which is diagonal in the flavor basis. These relaxation rates are similar in form to those 
proposed in refs. (281,150, 58, 59], within the context of active-sterile conversion or flavor conversion in supernovae. 

• The long time dynamics is approximately described by an effective Weisskopf-Wigner approximation with a 
non-hermitian Hamiltonian. The real part includes the "index of refraction" and the renormalization of the 
frequencies and the imaginary part is determined by the absorptive part of the second order self-energy and 
describes the relaxation. While this (WW) approximation describes mixing, oscillations and relaxation, it does 
not capture the dynamics of equilibration. 

• For time t >> rj~2 the two point function of the neutrino fields becomes time translational invariant reflecting 
the approach to equilibrium. The asymptotic long time limit of the one-body density matrix reveals that the 
density matrix is nearly diagonal in the basis of eigenstates of an effective Hamiltonian in the medium (|6.14[) with 
perturbatively small off-diagonal corrections in this basis. The diagonal components in this basis are determined 
by the distribution function of these eigenstates. 

• "Sterile" neutrinos: these results apply to the case in which only one of the flavored neutrinos is "active" 
but the other is "sterile" . Consider for example that the "muon neutrino" is sterile in the sense that it does 
not couple to the bath. This sterile degree of freedom is thus identified with the in medium eigenstate "2" 
because in the absence of mixing = its dynamics is completely free. The "sterility" condition corresponds to 
setting the matter potential = and the self-energy = with a concomitant change in the dispersion 
relations. All the results obtained above apply just the same, but with r^(fc) = 0, from which it follows 
that T2(k) = r ee (fc) sin 2 9 m (k). The final result is that "sterile" neutrinos do thermalize with the bath via 
"active-sterile" mixing. If the mixing angle in the medium is small, the equilibration time scale for the "sterile 
neutrino" is much larger than that for the "active" species, but equilibration is eventually achieved nonetheless. 
This result is a consequence of "active-sterile" oscillations which effectively induces an interaction of the sterile 
neutrino with the bathJH [H, US [H, H . 

Although the meson field theory studied here describes quite generally the main features of mixing, oscillations and 
relaxation of neutrinos, a detailed quantitative assessment of the relaxation rates and dispersion relations do require 
a full calculation in the standard model. Furthermore there are several aspects of neutrino physics that are distinctly 
associated with their spinorial nature and cannot be inferred from this model. While only the left handed component 
of neutrinos couple to the weak interactions, a (Dirac) mass term couples the left to the right handed component, 
and through this coupling the right handed component develops a dynamical evolution. Although the coupling to the 
right handed component is very small in the ultrarelativistic limit, it is conceivable that non-equilibrium dynamics 
may lead to a substantial right handed component during long time intervals. The study of this possibility would be 
of importance in the early Universe because the right handed component may thereby become an "active" one that 
may contribute to the total number of species in equilibrium in the thermal bath thus possibly affecting the expansion 
history of the Universe. 

Another important fermionic aspect is Pauli blocking which is relevant in the case in which neutrinos are degenerate, 
for example in supernovae. 

These aspects will be studied elsewhere. 
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